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We calculate 1-loop polarization in bilayer graphene in the 4-band approximation for arbitrary 
values of frequency, momentum and doping. At low and high energy our result reduces to the polar- 
ization functions calculated in the 2-band approximation and in the case of single-layer graphehe, 
respectively. The special cases of static screening and plasmon modes are analyzed. 



I. INTRODUCTION 

Graphene, a one-atom-thick layer of graphite, attracts 
a lot of attention of both theoreticians and experimen- 
talists since it's fabrication Quasiparticle excitations 
in graphene have a linear dispersion at low energies and 
are described by the massless Dirac equation in 2-t-l di- 
mensions. Theoretically such a behaviour was predicted 
long time ago and its numerous consequences were ex- 
perimentally checked after the discovery of graphene in 
laboratory. 

Contrary to the case of single-layer graphene (SLG) 
low energy excitations of the bilayer graphene (BLG) 
have parabolic spectrum, although, the chiral form of the 
effective 2-band Hamiltonian persists because the sub- 
lattice pseudospin is still a relevant degree of freedom. 
This changes many electronic properties of the material 
(for review, see 0, 3) compared to the case of mono- 
layer graphene. However, the low energy approxima- 
tion in bilayer graphene is valid only for small doping 
n < lO^^cm"^, while experimentally doping can attain 
10 times larrar densities. For such a large doping, the 4- 
band model @ should be used instead of the low energy 
effective 2-band model. 

In the literature, the screening effects in bilayer 
graphene were mainly^tudied within the low energy effec- 
tive 2-band model and in the presence of a magnetic 
field in [1, . Dynamical polarization plays an important 
role for finding plasmon excitations as well as for study- 
ing a gap equation and excitonic condensates both in 
single layer [l^ [ll| and bilayer graphene 0, . Some at- 
tempts to obtain analytical results in the 4-band model 
for bilayer graphene were performed in papers [l^ - [l^ . 
An exact calculation of the polarization function in the 
4-band model is interesting also from the pure theoret- 
ical viewpoint because then we can sec how the known 
results for the SLG [H [H] and the 2-band BLG [13 are 
recovered as limiting cases. 

Recently, a lot of attention is paid to investigate the 
properties of polarization operator in SLG |18l-|22|. The 
most general expression for dynamical polarization of 
SLG at finite temperature, chemical potential, constant 
impurity rate, quasiparticle gap, and magnetic field is 
given in Rcf. j23| . 



In this paper, we calculate the BLG dynamical polar- 
ization in the 4-band model within the random phase ap- 
proximation (RPA) for arbitrary wavevector, frequency 
and doping. Our results can be considered as an ex- 
tension of the results obtained in [l^l, although, those 
results were obtained in a slightly different approach. In 
Sec. |ll]we describe the model used and present our main 
result for the polarization function. We consider in Sec. 
nil Al thc static polarization function and compare it with 
the corresponding SLG and 2-band BLG results. In Sec. 
IIIIBI we focus on the long wavelength limit and study 
plasmons. Finally, we provide the details of our calcula- 
tions in Appendix A. 



II. RPA CALCULATION 

We model BLG in the Bcrnal stacking arrangement 
, where for two hexagonal lattices one sublattice of the 
bottom layer is a near-neighbor of the opposite sublattice 
of the top layer. In the tight-binding approximation, we 
have the following Hamiltonian: 



H 



(2.1) 



k a, 13 = 1 



Here Vfc = (al (k) , {k) , (k) , (k))^ , a°"'{k) and 
h°"^{k) are destruction operators of the Bloch states of the 
two triangular sublattices on the graphene layers a = 1 , 2 
with additional fiavour index a that encodes spin and 
valley. Further, is the electron density on layer a, 
Viiik) = V22{k) = 27re^/(Kfc) is the Coulomb interaction 
of electrons on the same layer, and electrons on different 
layers interact via Vi2{k) = V2i(fc) = Vii{k)e~^'^ , where 
d is the distance between the layers and k is dielectric 
permittivity of the substrate. One-particle Hamiltonian 
has the following form: 
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where t±^ ^ 0.4eV^ is the inter-layer hopping amplitude, 
e{k) = hvpikx +iky), and vector k = (k^, ky) describes a 
deviation from the if (^ = 1) and K' (.^ = — 1) point in 
the Brillouin zone [l^. Below we will consider only K- 
valley. The one-particle Hamiltonian can be diagonalized 
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with the help of unitary matrix U. Then one obtains the 
foUowing 4-band spectrum: 



H, = U^'dmg{E+, E^, -E^)Uk, 



Ef = J{hvFkY+tl/A±t^/2. 



(2.3) 
(2.4) 



In what fohows, we put for simphcity t± = t and rescale 
all momenta k k/hvf. Then the Fermi momentum 
equals kp = ^ //(/i + t) and charge density at zero tem- 



perature is given by n 
density n 

the higher density n = 10^'^ 



fc|,/t^lO"cm-^ The charge 



10^^ cm ^ corresponds to [iji = 0.1, while 



corresponds to kp = t 



and [ijt « 0.6. Here [i is chemical potential (Fermi 
energy). In the usual units, kp ^ t corresponds to 
kp « 0.061-1. 

If we denote the polarization matrix as follows 
-i{pa{io,k)pi3{-Cj,'k)) = S'^^^k - k)2Uai3ik), then the 
interaction potential in the RPA is given by 



V,fik) 
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(nii + ni2). 
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Here a = /{hvpn) is the effective coupling constant 
in graphene and d k, 3A is the distance between the 
graphene layers, which is relatively small, so we can 
set in all exponents d = (even for the largest possi- 
ble momentum e^^""^ « 0.85). Then we have = 
= l/ke{uj,k), where dielectric permittivity equals 
e{uj, k) ^ 1- aH(w, k)/k with H = 2(Hii 4- H12). The 
system is degenerate with respect to spin and valley de- 
gree of freedom therefore further we will consider polar- 
ization for one flavour degree of freedom: H -> H/7V/, 
Nf ~ 4. Then 1-loop polarization is given by: 



n— — oo 
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Summation over Matsubara frequency can be easily done 
[2?! . Then performing the analytic continuation through 
the replacement iujm — > w + iO, the retarded polarization 
function reads 
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where indices A and a denote bands and Fij is a 4 x 4 
matrix responsible for the chiral structure. It is defined 
as follows: 

F,,{q,p)=TT{Z-^A,ZA,), Z^U-^Up, (2.10) 

and Aj is diagonal matrix with all zero elements except 
unit at position j. We can find 



F{q,p) = 



where 



U++ V— V-+ u+- 

y— U++ U+- V-+ 

v+- U-+ u— V++ 
U-+ V+~ V++ u— 



(2.11) 



interband ones. At zero temperature, the Fermi func- 
tions in Eg. p.gp reduce to simple step functions. Then 
our retarded polarization can be presented in the follow- 
ing form: 

H(w,fc) =H°(cj,fc)+H+(w,fc)+H-(w,fc), (2.14) 
where 

H"(a.,fc) = / 5: ("^S^^^^C-^t/-- 



u^-[E[-UEt'? 



(2.15) 



1 + sw 
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e['^ = ^q^ + t2/4 + st/2, and 9qp is the angle between 
vectors p and q. Diagonal elements of F describe in- 
traband transitions while off-diagonal are responsible for 



(2.12) 



(2.13) 
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Clearly, H° does not depend on chemical potential and 
characterizes the polarization at zero doping. It gives 
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the main contribution to screening. The functions H"*" 
and n~ incorporate the effects of doping and are mainly 
responsible for plasmon modes. It is obvious that 11+ 
can be evaluated immediately if n~ is found for arbitrary 
values of fi and t. We have 11+ j = — t) n~_^ 

Let us comment on the chirality matrix p. lip . In the 
two limiting cases of weak {t — > 0) and strong (t — > oo) 
couplings when the spectrum reduces to Eg = q and 
Eg = q^/t, respectively, the chirality matrix F is strongly 
simplified and depends only on one parameter Uq. Then 
the polarization function per one flavour degree of free- 
dom [2^ equals 



(Pq 1 — Uq, 



q.k 



Eq + Eqj^k 



(Pq 1 + Uq^k 



2 a;2 _ 

Eq+k — Eq 



E. 



q+k) 



C^2 - {Eq - Eq + kf 



(2.17) 



where Uq^k = {'^T^HqHq+k) / {'^EqEqj^k) which is equal to 
cos 9q.q+k and cos 29q^q+k for weak and strong couplings, 
respectively. Note that at weak coupling 11+ = n~ while 
n+ = at strong coupling. 



In what follows we will consider an intermediate case 
for which fj, < t (only this regime is experimentally rele- 
vant). In this case, 11+ = 0. Calculation of II'^ and n~ 
is straightforward and the result can be written down in 
the following compact form: 



n{oj, k) ^ n°(a., k) + n- {lu, k) = - ^ - ^^^^^ + ^" - c^g. + 9t-. + 9t+. , (2.18) 



where 



with 



2lo 



3k^ - k^ + 5ij') + 2cj' 



2uj 



{u? -k^ -t^\-B [a;(cj + 2t) - k^\ ) 



2(a;2~fc2)^ 



;2 - fc2 



= Vw2 _ fc2 log ^^sgn(fc2 + y/k^ 



-, /i* = 2/i + t - cj, gu, = 



aj2 — fc2 0^2 — k^ 



\//c2 _ a;2 



tan" 



.1 Vfc2 - w2 



(^* ^ t - w) 



(2.20) 
(2.21) 



uj{t + a;)| log y 4- n 



fc2 ± ojr 
up — k"^ 



pI - r{uj + t) 



t)\ 



(2.22) 



Here expressions iO(...) are responsible for choosing the 
correct branch of the cuts. The square root and log- 
arithm have a branch cut discontinuity in the complex 
plane running from — oo to 0. Equations (j2.18p - (j2.22p 
are our main results. Details of the calculations as well 
as the expressions for the real and imaginary parts are 
given in Appendices lA II and lA 21 In the weak coupling 
limit t — >■ one can easily reproduce the results obtained 
in [l^ up to the overall factor 2 which reflects the bilayer 



nature of the system (in this case we should formally as- 
sume that fi > t and take into account n+(aj, k)). In the 
strong coupling limit t ^ fi,k,u!, in order to reproduce 
the results obtained in [13] one should take into account 
terms of order Ek = k'^ /t. 



It is convenient to normalize polarization with respect 
to the density of state at the Fermi level D{^) = Nf{t + 
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FIG. 1: Normalized polarization function for /i/t = 0.6. Panels (a) and (b) show density plots of the real and imaginary parts 
of the normalized polarization bubble defined in Eq. (|2.18p . respectively. Panels (c) and (d) present constant frequency cuts 
for u)/^ = 0.5, 1.0, 1.5, 2.0, 2.5. 



2//)/47r. So we introduce normalized polarization: 

n(.,fc).-2"°(-'^) + "'(-'^). (2.23) 
t + 2/i 



that the corresponding plots are very similar to those in 
Q . The static case a; = and the long wavelength limit 
— )■ are considered in Sees. IIII Al IIIIBI 



Finally, dielectric permittivity in terms of normalized po- m. ANALYSIS OF TWO PARTICULAR CASES 

larization is given by: 



e(w,fc) = K(l + 27raZ?(M)2^^ 1 , " " 



k 



A. Static screening 

The static limit a; — > is relevant for screening 



(2.24) of charged impurities. Performing some mathematical 
In Fig. [T]we plotted n(w, k) for fi/t = 0.6. One can note transformations we find that Eqs. ()2.18|) - ()2.22p imply 



n(uj = 0, fc = - log u ; tan^i - - Vt^ - k^ tanh"^ 



(3.1) 



VW^HKt^) f'^ + _ ^ COS- 4±^ ~ ttanh- + ^ ^t^^-Mt + M)] 



2k 2/x 



2k 



Vfc2 + 



fc2 - 2t^i 
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yW—k^ sinh 



fc2 




/c2 + t2 



The behavior of normahzed static polarization n(fc) = 
— 2II(a; = 0, k)/{t + 2u) and the corresponding polariza- 
tions for monolayer [l5| and bilayer in the 2-band ap- 
proximation [l^ are shown in Fig. [2^-[2t as functions of 
the normalized momentum k/kp. We see that the po- 
larization function calculated in the 4-band model has 
a discontinuity at k = 2kp similar to that found in the 
2-band model (see Fig. however, it does not go to 
a constant value at large momenta. Rather it grows lin- 
early as in the case of monolayer graphene (see Fig. [2f). 
For fi/t — >■ 0, the polarization function is similar to the 
polarization function in the 2-band model [l^l and tends 
to the SLG polarization function for /i/i ^ 1. Dielectric 
permittivity at large k for bilayer graphene in the 4-band 
model equals e{k) = l + naNf /A, whereas e(fc) = 1 in the 
2-band model. Note that e(fc) = l-|-7ra7V//8 for the SLG, 
therefore, we conclude that permittivities in the BLG in 
the 4-band model and the SLG coincide in view of the 
replacement Nf —J- 2Nf for the BLG due to doubling 
of the number of layers. Since the static polarization dc- 



(a) 



n(k) 



2.5 
2.0 
1.5 
1.0 
0.5 
0.0 



(c) 
3.0p 
2.5 
2.0 
n(k) 1.5 
1.0 
0.5 
0.0 L 



y:j/t=0.01 



M=o.i 



k/kr 
///t=0.6 







n(k) 



(e) 

2.5 r 

2.0 
1.5 . 
1.0 
0.5 

o.ot 



k/kf- 
-band BLG 




2 



FIG. 2: The static polarization. Panels (a), (b) and (c) show 
plots of the normalized static polarization given by (|3.1|) at 
fi/t = 0.01, 0.1 and 0.6, respectively. Dotted lines correspond 
to the asymptotic values n(fc) = 7rfc/2(t -I- 2/x). In panel (e) 
and (f) we show, respectively the static limit of the polariza- 
tion function for the bilayer graphene in the 2-band approxi- 
mation obtained in [ij and monolayer graphene calculated in 
the Dirac approximation in [l5| . [25l| . [3| . Dotted line at panel 
(e) corresponds to the asymptotic value Tl{k) = log 4 while 
the asymptotic at panel (f) is n(fc) = nk/S/j.. 

pends only on the absolute value of momentum, the RPA 



improved Coulomb potential is given by the following for- 
mula: 



Vir) = / dk 



I- 



kJo(kr) 



k + 2TTaD{^jL)Ii{k) 



(3.2) 



At finite doping the polarization function has a discon- 
tinuity at k = 2kp, therefore, at large distances the po- 
tential behaves as 

, ^ , , 1 sui(rkp) , ,„ „. 

V{r) V-^. rkp ->oo. (3.3) 

r rkp 



For zero doping, the discontinuity is absent and leading 
asymptotic is determined by the long-wavelength behav- 
ior of the polarization function. We find 



V(r) ^ ri — > oo . 



(3.4) 



The RPA improved Coulomb potential is shown in Fig. 
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FIG. 3: The RPA improved Coulomb potential at finite and 
zero doping. At panel (a) dashed line corresponds to the 
potential at fi — OAt and solid line to zero /i. Panels (b) and 
(c) show asymptotics of the RPA improved potential at zero 
and finite doping, respectively. 
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B. Plasmons 

The polarization function in the long wavelength limit 
A: ^ t is given by the following expression: 



^i^^k)^ — [^I + t + —log 



2fi + t 



2fj, + t + uj 



uj(u} + 2t) , 2n 
log 



uj(uj ~ 2t) , 2t 
log 



t 



A{t + uj) ° 2t + uj A{t-uj) 



2/x 



- LO 



t + UJ 

(3.5) 
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The plasmon dispersion relation is determined by the 
equation e{k,uj{k)) — which immediately gives: 



L0{k) 



2t 



(3.7) 



This is the general expression for the plasmon mode in 
2D systems, which for general spectrum of quasiparticle 
excitations can be written as [17| : 



w(fc) 



2k dq 

Equivalently this formula can be written as: 



Uj{k) = 2TTi 



(3.8) 



(3.9) 



where n = Nfkp/ATT is actual two dimensional density of 
particle while D{^) is density of states at Fermi level. In 
the case of the SLG D{^i) - so a;(fc) - fci/^ni/*. 

We solve equation Re[e(fc, cj(fc))] = numerically for 
free-standing graphene (i.e. k = I). Results arc shown at 
Fig. m We see that except " classical" plasmons with low 
energy behaviour (j3.8p we also have modes with linear 
behaviour and high- energy modes that are analogous to 
the TT-plasmons [26[. The corresponding dispersion rela- 
tion for small momenta are: 



uj{k) 
uj{k) 



2k{t + fi) k^t^^{t + fi) 



t + 2fi 

t + —^k\og 
2k 



it 



2^,Y 

2fi 
t 



1 + — 



(3.10) 
(3.11) 



However, contrary to the "classical" plasmons these 
modes cannot be considered as fully coherent collective 
modes, because they lie in the highly damped area which 
corresponds to the grey filling on the plots. Bound- 
aries of damped area are determined by the equation 
Imn(fc,w(fc)) = which can be easily solved, and we 
obtain 



w±(fc) = \lj 



{k±kFY - 



t 

2^^ 



(3.12) 
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FIG. 4: On the panels (a), (b), (c) we present dispersion rela- 
tions (black solid lines) for plasmons in free-standing graphene 
at densities ji/t = 0.6, yi/t = 0.3 and ji/t = 0.05 respectively. 
Black dashed lines describe classical plasmon (|3.8|l and high- 
energy plasmon (|3.1ip . Black dotted lines describe additional 
low-energy plasmon given by (|3.10|l . Filled areas show do- 
mains with nonzero imaginary part of the polarization whose 
boundaries (red dotted lines) are given by equations (|3.12p . 



that describes boundary of the single particle excitations 
continuum (Landau damping). Note that contrary to 
the normal 2D electron gas, plasmons damp at smaller 
momenta due to the interband transitions. 
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IV. CONCLUSION 

In this paper we have derived a compact analytic 
expression for the dynamical polarization for bilaycr 
graphene in the 4-band model in the random phase ap- 
proximation. Our results are valid for arbitrary values of 
wave vector, frequency, doping and interlayer coupling. 
Analysing the polarization as a function of the interlayer 
coupling we recovered the expressions for the monolayer 
graphene polarization (weak coupling) as well as for bi- 
layer graphene in the 2-band model (strong coupling). In 
the case where doping is smaller than the interlayer cou- 
pling we found the polarization function in the static and 
long- wavelength limits. Using these results, we have ob- 
tained the RPA improved Coulomb interaction and the 
dispersion relation for the plasmon mode. 

We put aside temperature effects and effects of the 
finite distance between layers, however within such for- 



malism they can be easily investigated and we postpone 
this investigation for a separate publication. 
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Appendix A: Calculation of the polarization function 

In this appendix we present some major steps in the calculation of the normalized polarization function. All 
quantities are evaluated in the units of energy (Sec. |ll]). We restrict our consideration to the case a; > because the 
polarization function for negative uj can be obtained through complex conjugation. 



1. n''(a;,fc) calculation 

In order to calculate n*'(a;, k) given by p.lSp . it is convenient to introduce the following variables: 

y = Eq+ Eq+k, Z = AEqEq+k ■ 

Then the measure of integration transforms as follows: 



7r/2 

j qdqde = 2jdq^Jde^^ J 



dy 



zdz 



(Al) 



(A2) 



{y-'-ef + ee 

y2 _ 



Performing integration over z, we get 

nO(w,fc) = 



2 \ w^-{y-t)^ w^-{t + yy 



(y^ — /c^)^^^ (w'^ — y"^) 



where (511° (w, k) is obtained by the proper change of the variables 
Jn°(a;,fc) = 



(A3) 



dy I V^^jy - t) ^ Vy^jt + y) y {3k^ - k^ {t^ + 5y') + 2y^) \ ^ ^ ^^^^ 



dy uP' + yt — k"^ 
2 w'^ — y^ 



(A5) 



y VW+k^-t -t/2-y/Wj4+k^ -t/2-y/W/4+k^ J 
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Now we can easily calculate the imaginary part for a; > 



Imn"(c.,fc) ^ / 3fc4 - P + 5^^) + 2w^ _ \e-oj{Lo-t)\ + \P-Lo{u + t)\ \ _ 7^^^)+ fA6) 



The real part is calculated treating all divergences in the principal value sense. After some algebra we obtain 



■ tan 



t 



(A7) 



4i:j 



log 



(/c2 + - wV)' 



4:UJ 



loe 



(fc2 + i2_a;2) (/c2 _^(2i + a;)) 



{k^ - uj{t + u;)Y 



+Re 



In the i — >■ limit, we get 



lim n°(w, k) = -- , — , 

t^o 4 Vfc2 - t^2 _ jO ' 

where we performed the shift a;2 — > ^2 + iO in order to reproduce the correct imaginary part. 

In the large t limit we must hold terms of order k'^ /t = Ej,. Then t appears only as an overall factor 



(A8) 



n°{uj,k) 

1/2 



loe 



El 






^El 


-4w2 





{Eu - w)' Eu + 2w 



{Eu + w)^ - 2a; 



iTT ( ( 1 - 1^") [a; - £;fc] - ('l - 1^ ) [2a; - £^^1 



(A9) 



2. n (cjjfc) calculation 



In order to calculate 11 (a;,fc) given by (j2.16p we introduce new variable r ^ Eq — t/2. Then performing some 
algebraic manipulations we find 



M 27r 



u-i.,k)^J'^J. 



IGtt V J' + r — Lo 



8r-4t , 



(AlO) 







where 



k^ - 2A:2 (2r2 + 2ra; + Ijj{uj - t)) + (2r + a;)2(t - a;)2 (fc2 - (2r + a;)(2r + t + Ijj))' 



fc2 + 2fc^r(r + i)cos + (2r + u;)(i - w) A:2 + 2kyjr{r + t) cos - a;(2r + < + 



(All) 



All divergences should be dealt with the prescription a; — )• w + iO. Then we can integrate over the angle using the 
following integral: 



1 

2tt J a + zeO + cosc 





sgn[a]6'(a2 - 1) sgn[e]6'(l - a^) 
— i- 



(A12) 
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We obtain the real and imaginary parts of the polarization function: 



\ 4(r + uj) 4(7' - Lu) 



-2r-t] , IniH" (w, fc) = - / — 



dr fRe[gi{uj)] Re[gi{-uj)] 



8r \ r + oj 



gniuj) = \/ (fc2 + (2r + uj){t ~ oj)f - Ak'^r^r + t)sgn {k^ + (2r + u)){t - w)) 



(A13) 
(A14) 



- (2r + w)(2r + t + w))%gn (fc^ - a;(2r + t + uj)) 
^ (P _ t^(2r + t + cj))^ - 4fc2r(r + t) 



(fc2 - (2r + a;)(2r + t + a;))%gn (r + § + w) 



= 4fc2r(r + t) - + (2r + w)(t - a;))^sgn ( r - | + cj ) + 

^ ^ y'4fc2r(r + f) - (P -tj(2r + < + w))^ 

(A15) 

We can calculate all integrals separately keeping regularization e of possible divergences at r = 0. In order to 
write down the answer in a compact form, we introduce the following notation. For any given function f{x), one can 



construct a new function f{x) 



by the following rule [29j: 



f{x) = sgn(fe - x){f{b) - f{r)) - sgn(a - x){f{a) - /{x)) 

a 

Then one can present the polarization in the following form: 



(A16) 



n, {uj,k) = -/i--iog — + ITT — e{fi~uj){0{pt^j~0{~pj) , (A17) 



where = Rf. + = + lui, and 



2uj 



LO — t 



Jt-u 



t-u, 2(^2 - 



2fi—uj 



J 



t+UJ 



LU — t 



(A18) 



+ {{lj - ty - k')ut-^ , 



: (A19) 



(A20) 



Here 

Puj = \/ fc^ 

and 



P-t2 



t — UJ 



2(^2 - P) 



2(1— UJ 



uZ(yj) - {{uj -t) - k )ut-uj{t - uj) 



2^-\-uj 



(A21) 



, /i = Ai + t/2, /,^(r)=tan-i 



pi sin (r/pi^) 



A't(r)-sgn(c.(w^-P-i^))/^(r) 

(A22) 



X\/pI-x- 



4Vk^ - uj^ 



, ^^^(a;) = sgn(a;(fc^ - a;^))t;^(x), u^{x) = ^^-j=^l£^, u^{x) ^ sgn{uj{k'^ - uj'^))u^{x) . 



4Vk^ - ^2 



(A23) 
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Expression (IA17P should be understood in the hmit e — 0. Taking this hmit expUcitly, we find 

n (w, k)^-n+ + I (A24) 



where 



prcg _ .\k'^ + (t — U!)oj\ / ~2/i-w,-fc^/w ^2/i+tj,fc^/Lj ^t+Lj.fc^/ti; ^2;i+w,fe^/(tj-t) ^2^i+uj ,k'^ / (ui-t) ,k'^ / (ui-t) 



(A25) 



^ hsgn(-)/.„. j - 2 1087^ pS-ft-^)^ 



rreg 1^ w)ti;| /„2/i+Lj,— w ^t+w,— ^2/i+i»;,t— w . /^uj,t—bj ^2/i— w.w . ^2/i+tj,t— / A 0«^ 

|/c^ + (i — a;)^! / 2 2\ -1 f ^ ^ ^\ ^ 



2^ ^0(fc^-c.^)cos'^ - j - - /r_ J^) + ^ t - c.) 

and 

Co = sgn(6 - a) (/- (6) - A'^ (a)) , C,^,, = sgn(6 - a) (/^ (5) - (a)) . (A27) 
In the weak couphng hmit < — >■ 0, we have the fohowing expression: 

n {uj, k) = -/i H , (A28) 



nt=0 _ 

ii,,, — 



n ((fc2 - t^2^ (fc2 _ ,^2 _^ 2w/i)) ((2^ - a;)yp~~(2^7^^ + P (sin"! (fc/w) - sin'^ ((w - 2^)/fc))) 

4VA:2 - u;2 

(A29) 



^ (2/.-c.)Vfc2-(2M-c.)2 fc2 (sin-i((,^-2M)/fc)+sin-i(^/fc)) 

=^g^(^--)(^ i7F=^f i7F=^f ij- (^30) 

We can unite real and imaginary part in one expression: 



" (^'^) = '^' 2(c.2-fc2) + 4 ' 



where 



" ^"+* 2(a;2-)t2)2 + 'Y\IJ^^ + C.2 - fc2 + 

(A32) 

k^-ujjt + cj) plp^\{uj'^~k^){uj{uj + 2t)-k^)\ iTr\k^ -ujit + u)\ ^ 
+ 2^ (fc2-c.(^ + t))4 + -k -t]-e[ui. + 2t)^k]), 

and the functions Guj and Qlji u ^^'^ determined in Eqs. (|2.20l) - (|2.22l) . 



[1] K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, Y. Zhang, S. V. Dubonos, I. V. Grigorieva, and A. A. Firsov, 
Science 306, 666 (2004). 



11 

[2] G. W. Semenoff, Phys. Rev. Lett. 53 2449 (1984); P. R. Wallace, Phys. Rev. 71, 622 (1947). 

[3] V. N. Kotov, B. Uchoa, V. M. Pereira, A. H. Castr o Neto and F. Gu inea, „arXiv:1012. 3484,. 

[4] S. Das Sarma, S. Adam, E.H. Hwang and E. Rossi, larXiv:1003.473TV 2. 

[5] Ed.McCann and V.I.Fal'ko, Phys. Rev. Lett. 96, 086805 (2006). 

[6] E. H. Hwang and S. Das Sarma Phys. Rev. Lett. 101, 156802 (2008). 

[7] R. Nandkishore and L. Levitov. Phys. Rev. Lett. 104, 156803 (2010) 

[8] E.V. Gorbar, V.P. Gusynin, and V.A. Miransky, JETP Letters, 91, No. 6, 314-318 (2010). 

[9] E. V. Gorbar, V. P. Gusynin and V. A. Miransky, Phys. Rev. B 81, 155451 (2010). 

[10] O.V. Gama yun, E.V. Gorbar and V.P. Gusynin, Phys. Rev. B 81, 075429 (2010), 

[11] J. Gonzalez larXiv: 1103.3650] 

[12] G. Borghi, M. Polini, R. Asgari, and A. H. MacDonald, 2009, Solid State Commun. 149, 1117. 

[13] G. Borghi, M. Polini, R. Asgari, and A. H. MacDonald, 2009b, Phys. Rev. B 80, 241402. 

[14] S. V. Kusminskiy, J. Nilsson, D. K. Campbell, and A. H. Castro Neto, 2009, Europhys. Lett. 85, 58005. 

[15] B. Wunsch, T. Stauber, F. Sols and F. Guinea, New. J.Phys. 8, 318 (2006). 

[16] E. H. Hwang and S. Das Sarma, Phys. Rev. B 75, 205418 (2007). 

[17] Sensarma, R., E. H. Hwang, and S. D. Sarma, Phys. Rev. B 82, 195428 (2010). 

[18] P. K. Pyatkovskiy J. Phys.: Condens. Matter 21, 025506 (2009). 

[19] A. Qaiumzadeh and R. Asgari, Phys. Rev. B 79, 075414 (2009). 

[20] K. Shizuya, Phys. Rev. B 75, 245417 (2007). 

[21] R. Roldan, J.-N. Fuchs, and M. O. Goerbig, Phys. Rev. B 80, 085408 (2009). 

[22] R. Roldan, J.-N. Fuchs, and M. O. Goerbig, Semicond. Sci. Technol. 25, 034005 (2010). 

[23] P. K. Pyatkovskiy and V. P. Gusynin, Phys. Rev. B 83, 075422 (2011). 

[24] V.P. Gusynin, S.G. Sharapov and J.P. Carbotte, Int. J. Mod. Phys. B 21 4611 (2007). 

[25] E. V. Gorbar, V. P. Gusynin, V. A. Miransky and I. A. Shovkovy, Phys. Rev. B 66, 045108 (2002). 

[26] S. Yuan, R. Roldn and M. L Katsnelson. iarXiV;i l03.5350 

[27] Really T^^^__^(in„ + a) ^i(in„ - n + b)'^ = {npia) - nF{b))/{a - b) where npix) = {1 + ex.p{{x - fi) /T))-\ 

[28] In the weak coupling regime, the number of flavours effectively doubles. 

lb 

[29] One can easily note that f{x)\ = f drf {r)sgn(r — x). 



